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We study the propagation and transmission of Alfven waves in the context of cylin- 
(JZ ■ drical geometry. This approximates the polar cap region of aligned pulsar with strong 

iy~j ■ magnetic fields. Non-propagating region appears in the presence of rotation. The dis- 

placement current further prevents the low frequency modes from propagating near the 
,__! , stellar surface. The transmission rates to the exterior through the surface are calculated. 

The rates increase with the frequency and the magnetic field strength. The transmission 
r- — | also depends on the helicity states of the waves, but the difference becomes small in 

the high frequency regime. We also point out the possibility of the spin-up by outgoing 
wave emission in the low frequency regime, if a certain condition holds. 

Subject headings: Star: Magnetic Fields — Star: Neutron — Star: Magnetar — Gamma 
Rays: Bursts 
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1. Introduction 



* 



Soft Gamma Repeaters (SGRs) and Anomalous X-ray Pulsars (AXPs) belong to a rare class 
among thousands of neutron stars. Both objects are different from so far known neutron stars 
in their radiation spectrum and spin periods. SGR-like outbursts were also discovered in AXPs 
(Gavriil, Kaspi & Woods 2002; Kaspi et al. 2003), and hence there are similarities between two 
peculiar neutron stars. They are likely to be young and isolated, but have intense magnetic fields 
in the 10 14 - 10 15 G range. See Mereghetti & Stella (1995); Kouveliotou et al. (1998) for their initial 
observations and also Woods (2003) for recent reviews. Further evidence of the strong magnetic 
fields has recently come from spectral line feature(Ibrahim et al. 2003). These objects are called as 
magnetars (Thompson & Duncan 1993, 1995, 1996). The magnetar model has been able to explain 
the peculiar observational properties. See e.g. Thompson (2000); Lyutikov et al. (2002) as recent 
reviews. The burst emission and non-thermal X-ray radiation are supplied by the decay of strong 
magnetic fields. Twisted magnetic fields relevant to the activities are transported from the core to 
the surface by ambipolar diffusion in the magnetars model (Thompson et al. 2002). 

The Alfven waves are likely to be excited by sudden disturbances in various magnetized objects. 
Short duration bursts in SGRs may be associated with starquakes driven by magnetic stress, since 



the similarity in statistical property is reported between the SGRs and earthquakes (Cheng et 
al. 1996; Gogus et al. 1999, 2000). The Alfven waves are very interesting in their property and 
physical mechanism. The fluid displacement is perpendicular to both of the wave vector and the 
magnetic field, and restoring force arises from the magnetic tension of the field lines. The plane 
Alfven waves of finite amplitude can propagate at constant speed in a homogeneous incompressible 
medium without any distortion of the waveform. This contrasts with sound waves, which may 
steepen to form shocks because of the non-linearity. As for the physical mechanism, torsional 
Alfven waves can transport angular momentum, and account for the spindown of rotating objects, 
such as magnetized stars with convective envelopes, interstellar clouds threaded by the Galactic 
magnetic fields and so on. See e.g. Shu (1992) for details. 

We consider the Alfven waves produced by shaking the magnetic field lines in a rotating 
magnetar in this paper. The origin of abrupt disturbances relevant to the bursts is not addressed 
here, but the subsequent propagation and ejection to the exterior are examined. The total energy 
in shorter duration bursts is AE < 10 41 ergs, and is a small fraction of available magnetic energy 
~ 10~ 5 Eb, where the total magnetic energy Eb is estimated by the expected dipole field strength 
B = 10 14 — 10 15 G. On the other hand, giant flares with AE ~ 10 44 ergs should involve more 
drastic change of magnetic field configuration on a global scale (Ioka 2001; Woods et al. 2001). 
We here use the linearized perturbation equations, which are applicable to less energetic events 
AE /Eb oc (SB/B) 2 <C 1, in the shorter duration bursts. Our results of the wave propagation may 
not be applied to the giant flares, but there may be some similarities in nature even for 5B ~ B. 
The ejection of Alfven waves to the exterior was already estimated so far, but our model is beyond 
a simple model. We here explicitly calculate the transmission in a cylindrical geometry. We also 
calculate the energy and angular momentum extracted by the waves through the stellar surface. In 
section 2, we provide a geometrically simplified model. The propagation in the neutron-star crust 
and transmission to the exterior can be examined in a concrete way owing to the simplified model. 
The shear in the solid crust acts as additional restoring force. Therefore, the shear-Alfven waves 
may be more adequate for the name. The transmission rate is numerically calculated for the waves 
with higher frequency. We discuss the implication of our results in section 3. 



2. Alfven Waves in Rotating Cylinder 

2.1. Model 

In this paper, we study the propagation of the shear-Alfven waves through a neutron-star 
crust and the transmission to the exterior. The star is assumed to rotate around z-axis with a 
constant angular velocity £1. The magnetic field is uniform along the z-axis, i.e. B = Bqs z , where 
Bq is a constant. Our consideration is limited to a cylindrically symmetric slab region within the 
radius w*. The magnetic fields are rather easily incorporated in the cylindrical model. Carroll et 
al. (1986) studied the oscillation spectra of the magnetized/unmagnetized cylindrical stars. Their 



results show that the periods of the torsional and shear modes are in good agreement with those of 
spherically symmetric unmagnetized star. This suggests that the approximation is good for some 
modes. 

The cylindrical approximation may be adequate for the polar cap region and the interior for 
an aligned rotator. In this case, the radius m* is given by the size of the polar cap, zu^ < R sin 8 P = 
a/-R 3 S1/c ps 10 (T/ls) -1 ' 2 cm, where R ~ 10 6 cm is stellar radius and T is the rotational period. 
The exterior magnetic field lines originated from this region are extended to infinity. On the other 
hand, the field lines outside the polar cap become closed ones. The plasma along the closed lines 
corotates with the central star. In our cylindrical model, we consider the shaking of open field lines 
only, which may lead to particle acceleration, magnetic reconnection, wind, and eventually to the 
radiation. Explicit treatment of the energy transfer to the radiation is also beyond the scope of 
this paper. 

The shear stress in the crust also acts as a restoring force, and is coupled to the Alfven waves. 
The solid crust ranges from z = —d to the surface z = 0. The depth d of the outer crust is given by 
d ~ 10 5 cm. Our consideration is limited to thin outer crust only. The slab geometry is valid there 
because d <C R. Gravitational acceleration in this region is almost constant and is approximated as 
—ge z , with g = 10 14 cm s -2 . The density distribution p(z) is given by the integration of hydrostatic 
equation for neutron star crust (\z\ < d), where degenerate electron pressure is dominated. The 
analytic expression is available as (Blaes et al. 1989) 

p = 8.0 x 10 1 [0|/cm) + 2.5 x 10" 4 (|z|/cm) 2 ] 3/2 g cm" 3 . (1) 

This explicit form will be used for the numerical calculations in the subsequent sections. 



2.2. Wave equation 

We consider the propagation of shear- Alfven waves based on the linearized perturbation theory. 
In the deep interior, gravitational force is so large that the displacements to the vertical direction 
are not easily to be induced. The horizontal displacements of the disturbances are likely to be 
dominated, and they are coupled with the Alfven waves. Some restrictions are imposed on the 
displacement vector £ of the materials in order to extract the waves coupled with shear and magnetic 
stress. The waves propagating along the z-axis are assumed to satisfy £ z = V • £ = 0. This means 
that the waves are transverse and are decoupled from the compressional modes. The Lagrangian 
perturbations of density and pressure are zero because they are proportional to V • £, i.e. Ap/p = 

-v • i= o, Ap/ P = -{diogp/diog P ) ad v • i= 0. 

The general forms of the displacements can be expressed by the Fourier and Bessel functions 
with respect to time and cylindrical radius, respectively. We further simplify the displacements by 
assuming nodeless functions in cylindrically radial direction, i.e. neglecting the radial structure. 



The regular form near the z-axis is simply given in the cylindrical coordinate (zu, (f>, z) as 

£ oc (e ro ± te^)m m - l ti± m M)e~^ tTm(f,) , (2) 

where m is a positive integer, m > 1. The mode for the upper sign in eq.(2) is said to have 
positive helicity, whereas that for lower sign has negative helicity. We may limit the Fourier mode 
to a certain frequency region, using a general relation £± m _ w = £* m u , where * means complex 
conjugate. We define the frequency in corotating frame a± = to =F m Q f° r each helicity state £± 
and limit the range to a± > 0. From now on, we will omit the suffixes mw in order to avoid 
complicated notations. The displacement (2) satisfies (V x £) 2 = and is therefore decoupled 
from the vorticity as it is desired. The displacement corresponding to m = 1 has clear meaning, 
and can be written in the Cartesian coordinates as £ = (e x ± ie y )t;±e~ lu>t . This mode represents 
uniform motion in horizontal direction. A slightly different treatment is necessary for the axially 
symmetric perturbation m = 0. The regular displacement satisfying above conditions is given as 
£ = ^ ( f,(t,w)zue c f > , which corresponds to the velocity perturbation 5v = (d^/dt)we^ = dQzue^. 
This kind of perturbation merely represents impulsive jump of angular velocity as O — ► Q + 8Q, 
and is no longer considered in this paper. 

The restoring forces for the modes are elastic shear stress 6 Si and electromagnetic force 5Fi. 
The linearized equation of motion can be written as 

O.Svi + (v ■ V)Svi + (Sv ■ V)v^\ = 5 Si + 5F U (3) 

where v = w^le^, and the relation between the displacement and velocity perturbation is 

6vi = ^i + (v-V)Zi-(g-V)vi. (4) 

The shear stress 5Si associated with deformation is 



5 Si = V, 



dxi dx l 



(5) 



with shear modulus n = 4.8 x 10 27 (p/10 11 g cm 3 ) erg cm 3 (Baym, &: Pines 1971). The 
electromagnetic force is given by 

SF = 5p e E + Pe 5E + -(5jxB+jx 5B). (6) 

c 

The crust is a perfect conductor, so that the electric fields for both unperturbed and perturbed 
states are induced by the material motion as 

E = — (vxB)= °-e w , (7) 

c V / c 



SE = — UvxB + vxSb). (8) 

The electric charge and current for the unperturbed state are 

Pe = ~- — ftB , j = p e v. (9) 

Z7TC 

These quantities (7) and (9) satisfy the force balance p e E + j x B/c = 0. These expressions are 
valid in the region smaller than light cylinder c/Q, or actual stellar radius R. Otherwise, we would 
need better approximation beyond the uniformly rotating cylindrical model. Therefore, our present 
model may be adequate to the behavior near the z-axis. As will be confirmed by the explicit forms 
of E and SE, the first term 5p e E(= (V • 8E/4n)E) in eq.(6) is proportional to (Qw/c) 2 and is 
smaller than other terms near the z-axis, i.e. for the region w 2 < £/£". We therefore neglect the 
first term 5p e E in the propagation of Alfven waves. Eliminating 5E by eq.(8), the Lorentz force 
(6) is reduced to 

SF = -(Sj-p e Sv)xB. (10) 

c 

The perturbation of electric current is determined by the Maxwell's equations: 

Vx5E = -~5B, (11) 

c at 

4ti- _ l d ^ 
V x5B = —Sj + -—SE. (12) 

c cot 

From eqs.(8) and (11), the perturbations of electro-magnetic fields are expressed in terms of the 
displacement vector as 

SB = B V Z £= Bom™- 1 ^-^ ± ie^e^^, (13) 



5E = ^ 



m—l r (\p\ l 

■ T^Tm^ile^i^ + fi^e-^^. (14) 

The stellar rotation induces the ^-component of the perturbed electric field, which is important for 
angular momentum transfer as discussed in section 2.6. Using these expressions and eliminating Sj 
by eq.(12), eq.(3) is eventually reduced to 



where 



47TC 2 p + £ 2 



2.3. Propagation 

The propagation of Alfven waves is studied in different contents. In a certain limit of pa- 
rameters, our basic equation (15) should be reduced to one previously studied. In the limit of 
non-rotating case, eq.(15) is reduced to one found by Blaes et al. (1989). The term V± is indepen- 
dent of the helicity state and is positive definite in this case. Another interesting limit of eq.(15) 
is obtained when both shear and relativistic effects are neglected. This limit corresponds to // — ► 
and c 2 — ► 00. In this case, the term V± becomes negative for a certain frequency region. Hence, 
the modes are non-propagating (evanescent). See e.g. Chandrasekhar (1961). We will discuss the 
evanescent property below. 

We consider the short-wavelength limit in order to see whether or not the waves propagate. 
The local dispersion relation is obtained by setting £± oc e in eq.(15) 

-A 2 + (wTmfi)(wTK-2)n) = 0, (18) 

where v is the velocity defined as 



( 47^ + B 2 
\4irpc 2 + Bl 



1/2 



(19) 



In the limit of \x — ► and c 2 —> 00, the classical Alfven wave velocity is recovered, i.e. v = 
Bo/^Airp. The shear acts as the restoring force deep in the crust unless Bq 3> 2 x 10 14 G. The 
displacement current contributes to additional inertia and is always important near the surface, 
where the propagation velocity becomes c. Solving eq.(18) for the frequencies a± = u> =F TnCt > 0, 
we have 



0± 



m 



2 



m 



2 



2 



i + ^(i_^) n + \ i + -(i-h) n 2 + v 2 k 2 \ . (20) 



1/2 



We do not have to consider the modes with negative frequency a± < because of the correspondence 
(£±j a ) ^ (£, —cr) as mentioned in section 2.2. The difference between two helicity modes is clear 
in weak magnetic field limit, i.e. vk/Vt — ► and h — ► 1. In this limit, the frequencies in the 



corotating frame are given by a_ ~ 2Q and a + ~ v 2 k 2 /(2Q). The mode £_ with frequency cr_ 
represents an inertial wave due to Coriolis force. On the other hand, the mode £+ with smaller 
frequency o+ (^ <r_) represents a drift wave. Our problem is limited to cylindrical geometry, but 
the similar oscillations are also possible in spherical geometry. See e.g. Levin & D'Angelo (2004) for 
comparison in the non-relativistic limit, although careful treatment is necessary to convert results 
between different geometry. 

Like the Alfven waves in classical treatment, the waves with low frequencies become non- 
propagating due to the Coriolis force. The mode £+ is always propagating as far as a + > 0, 
whereas the mode £_ is non-propagating for the frequency region < er_/f2 < 2 + m(l — h). 
The condition becomes < a^/O, < 2, when the relativistic effects are neglected. The evanescent 
property in our problem depends on the spatial position through h(z), which is calculated for the 
density distribution (1) and constant magnetic field strength Bq. It is general that the function h 
decreases from h ~ 1 in deep interior to h ~ near the surface with low density. The evanescent 
region therefore prevails near the surface. In Fig.l, we demonstrate the location of evanescent 
region using the local dispersion relation. The critical frequency a-/£l = 2 + m{\ — h) is shown. 
The depth of the evanescent region increases with azimuthal wave number m as well as the magnetic 
field strength Bq. 
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Fig. 1. — The critical angular frequency <j-/Vt = 2 + to(1 — h) as a function of the depth from the surface. 
Evanescent region for negative helicity mode £_ is given by the condition < <7-/Q < 2 + m(l — h). The 
critical frequencies are shown for m = 1,3,5 with Bo = 10 13 G and m = 5 with Bo = 10 15 G. Note that 
the evanescent condition becomes < cr^/Vl < 2 denoted by a dash-dotted line, if all relativistic effects 
are neglected. The evanescent region near the surface enlarges with increase of m and the magnetic field 
strength. 



The phase and group velocities of the shear- Alfven waves in the inertial frame are respectively 
riven by v p = uj/k and v g = dco/dk. The following relation is easily calculated 



to 



<J± ± 7TJ.S7 



w±i[2-m(l + /i)]n a±±±[2 + m(l-h)]n 



(21) 



This shows an interesting property that the phase and group velocities of the mode £_ are opposite 
in the propagating direction for the frequency range [2 + m(l — h)]£l < cr_ < mCl, which is possible 
for m > 3. The lower bound comes from the condition for the wave propagation. 

In Fig. 2, we show 1/N 2 = u; 2 /(c 2 /c 2 ) as a function of cr_/Q = uj/fl + m, using the local 
dispersion relation (18). The frequency corresponding to negative value of 1/iV 2 mean evanescent. 
The non-propagating frequency is located in low frequency region. We now consider the propagation 
of the negative helicity mode with 2 < a-/Q < m. The mode propagates as a wave in deep 
interior, and may leak out through to the exterior. The amplitude of the outgoing wave becomes 
small through the evanescent region. It is important to know the magnitude. The damping rate 
depends on the frequency and the depth d*. The exponential damping factor e~ K * is estimated as 
nd* « \N\wd*/c, which is not so large as roughly estimated ~ 10~ 4 for large m. This means that 
the damping of the amplitude is not so significant, and therefore the disturbances of some modes 
are penetrated into the surface. 
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Fig. 2. — Dispersion relation of negative helicity mode £_ as a function of angular frequency er_/f2 at 
z = 10cm from the surface. 



2.4. Exterior 

There are at least two possibilities as for the exterior of the stars. One is that surrounding 
plasma outside the star corotates with the same angular velocity. When the plasma is frozen 
in closed magnetic field lines, such a corotation is realized. However, the plasma freely moves to 
infinity along open field lines, so that the angular velocity is arbitrary in this case. We also consider 
this situation by assuming that the plasma exterior to the star is static in the inertial frame. The 
difference of the plasma motion in the background state is not so important for the high frequency 
modes of the perturbations (uj 3> Q), but crucial for the low frequency modes, as discussed below. 

When the plasma corotates with the same angular velocity fi, the perturbation equation cor- 
responds to the limit of eq.(15) with p = fi = 0. The wave solution is possible for the modes with 
(uj =F mO)(w ± 2S7) > and outgoing wave is expressed as 

£± = A±exp[ik±z], (22) 

where 

c 

We assume that the interior and exterior solutions are continuously matched, so that the boundary 
condition of the displacement at z = can be expressed as 

§ =tk±i± . (24) 

We next consider the non-rotating plasma in the exterior of the star. In this case, outgoing 
boundary condition at the surface is replaced as 

— — = i-£±. (25) 

dz c 

The shaking of the magnetic field lines induces electromagnetic radiation in the exterior vac- 
uum. Since 5B Z = everywhere, the perturbation should be matched with the transverse magnetic 
mode with the cylindrical wave guide (Jackson 1975). The perpendicular components of elec- 
tric fields, i.e. SE-^^SE^ are continuous. They are expressed by £ in the interior and vary as 
oc exp{iuj(z/c — t)} in the exterior. Thus we have the boundary condition (25) at the surface. 



2.5. Transmission coefficient 

In this section, we consider how much Alfven waves excited in the crust are transmitted to 
the exterior. We numerically solve eq.(15) from the deep interior z ~ 10 5 cm, corresponding to 



10 



the neutron drip in the density, to the surface. The following asymptotic solution is obtained in 
the deep crust (Blaes et al. 1989). It is a sum of outgoing and reflecting waves, which respectively 
propagate to +z and — z directions. We normalize the amplitude of outgoing wave to unity for 
simplicity and denote that of reflecting wave by R±. 



£± ~ M 4 { ex P 



— dz 



+ R± exp 



— dz 

v 



(26) 



Here, T± is an effective frequency of each wave mode defined as 

r ± = {(uTmfi)[wTK-2)fi]} 1/2 . (27) 

The outgoing wave solution in the exterior is given by eq.(22) with a transmission amplitude A±. We 
here assume that the exterior plasma corotates with Q. The interior solution should be continuously 
matched with the exterior one at the surface, z = 0. Thus, the transmission coefficients T± for each 
helicity mode are given by the relation 

T± = \A±\ 2 = 1 - \R±\ 2 . (28) 

Our numerical calculation is limited to m = 1. The lower bound of the frequency for the 
positive helicity mode is <r+ > 0, that is uj > O. On the other hand, that for negative helicity mode 
is limited by propagation condition cr_ > 30, that is uj > 2£l. In Fig. 3, transmission coefficients 
T± for each helicity mode are shown as a function of angular frequency in inertial frame to. The 
rotational effect may be important in the millisecond magnetar model. In the calculation, the star 
is assumed to rotate with a constant angular velocity, £l/(2n) = 10 3 Hz. This value is somewhat 
too rapid for SGRs, but we can regard it as an extreme one, which may be present at the newly 
borne phase. The rotational effect becomes small for much smaller value of 0,. For a comparison, 
we also show the result To for the non-rotating case ($7 = 0) by solid curves in Fig. 3. The difference 
between positive and negative helicity states is clear near the low frequency limit, i.e. uj < a few 
xQ. The difference however becomes almost negligible in the high frequency regime, say, uj > 10S7. 
The transmission rates gradually approach unity with increase of the frequency. The rates also 
increase with magnetic field strength, but T± ^ 1 in the low frequency regime even for 10 15 G. 
A significant fraction of waves suffers from the bounce at the surface, because of T± ~ 0.1 for 
uj < 10O. This means that the transmission time to the exterior is roughly estimated not by the 
crossing time d/v, but by d/v x T± . 



2.6. Energy and angular momentum loss 

We will calculate the rate of energy and angular momentum carried by the Alfven waves. 
The terms linear to the perturbation quantities vanish in the time average. Meaningful time- 
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Fig. 3. — Transmission rates as a function of angular frequency for the magnetic field strength 10 13 and 
10 15 G. The curves with labels T + and T_ denote results for the positive and negative helicity states. The 
solid curves with the label T denote results for the non-rotating star. 



averaged values come from square of the perturbation quantities. The time-averaged Poynting flux 
is therefore given by 



S = — 5ft ( SE x SB* 

8tt 



(29) 



where * means complex conjugate and ft() means taking the real part. We have here used a useful 
technique in averaging the product of complex quantities with the same harmonics time-dependence 
as shown in the textbooks, e.g. Jackson (1975). The power P radiated across the surface at z = 
is calculated by integrating S z over a circle within w* as 



P 



c 



2tt 



ft {dE^SB^ - 6E%6B„) zadwc 



(u T mtyB^ml 
Am 



2in 



ft -i£ 



: dz 



(30) 



z=0 



This is roughly estimated as P ~ (5B) 2 (irza^)c. The power Pr R \ ejecting to the corotating ambient 
plasma is evaluated using the boundary condition (24). The power Pr<n to the static plasma is also 
calculated with condition (25). The explicit expressions are 



(R) 



(u;Tmn)k ± Bfa 2 ™\Z±\ 
Am 



(31) 



-12- 

p ( 5 ) = i^ • (32) 

In a similar way, the angular momentum flux per unit time across z = can be calculated. 
The transport of the angular momentum results in the torque. The time-averaged torque due to 
the electromagnetic stresses can be expressed in terms of the surface integral. From the symmetry, 
non-vanishing angular momentum flux is z-component only, which is expressed as 



n z = -—I I a {6b;sb z + se;se z ) w 2 dw 



„, , , , » -i^±^ • (33) 

8(m+l)c 2 V dzj z=0 

The angular momentum transfer in eq.(33) comes from the electric part only, since SB Z = 0. The 
value is estimated as N z ~ (Qw^/c) x ((5S) 2 tu*(7rc[7 2 ), and is smaller by an extra factor (Qw*/c) 
compared with the rough estimate by the magnetic part only. Note that the angular momentum 
transfer comes from the magnetic term, i.e. B^B Z for non-relativistic case such as interstellar 
clouds. Depending on the boundary condition (24) or (25), the angular momentum flux is written 
as 

**<*> = 8(m + l)c 2 ' (34) 



(w T m n)ujnB 2 m 2 . m +2 \Z±\ 2 
8(m + l)c 3 



7V 2(s) o/„, i TV5 • y 6b ) 



It should be noted that the angular momentum flux can be negative by the outgoing waves 
into the static ambient plasma as shown in eq.(35). This is possible if the condition [to =F mQ)uj = 
a±to < is satisfied. That is, frequencies, a and to, are opposite in sign. This condition is the 
same as the criterion of gravitational radiation reaction instability (Friedman &; Schutz 1978). 
The oscillation mode counter-rotates when viewed in rotating frame with star, but corotates when 
viewed in the inertial frame. This provides a mechanism for converting the stellar rotational energy 
into gravitational radiation. The modes satisfying the condition {io^mQ)io < are very interesting. 
They carry negative angular momentum. The emission leads to the spin-up of the star. 



2.7. Order of magnitude 

There is a simple relation between the energy flux (30) and angular momentum flux (33) 
exerted by the Alfven waves: 

N z = m ^ W \ , R (36) 

2(m+l)c 2 V ; 
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This holds irrespective of the boundary condition (24) or (25). The total angular momentum 
loss/gain is related with the total energy radiated by the waves. During the waves emission, the 
stellar structure may change due to re-arrangement of magnetic fields, and the inertial moment I 
may be modified. For giant bursts in SGRs, such a possibility is discussed (Ioka 2001; Woods et al. 
2001). We do not consider the possibility here. The duration time of the wave emission is denoted 
by r. The change in the angular velocity AJ7 is expressed by the total energy AEa = Pt 

AQ_N z t_ mwl 20 / ™* \ 2 ( ^E A \ 

Q ' m 2(m + l)c 2 / A vio 4 c m ; Vio 38 erg S y' l ; 

where the total energy by the Alfven waves is not known, but is used for the typical burst energy in 
observed in SGRs. There are reports of glitches observed in AXPs(Kaspi et al. 2003), A17/J7 ~ 10 -6 . 
Our estimate is too small. Much larger total energy may be expected associated with the giant 
bursts. Even taking comparable energy observed in X-/7~rays in such events AEa ~ 10 44 ergs, the 
change is Ail /ft ~ 10 -14 , whereas |AJ1/0| ~ 10~ 4 in the giant burst. One reason of such smallness 
may come from too much idealized situation, in particular, our model giving 5B Z = exactly. This 
introduces a small factor (£lw*/c) ~ 10 -6 in the angular momentum loss N z . Moreover, larger size 
w* should be required also. Otherwise, a wind of particles and MHD waves or interaction between 
the crustal neutron superfluid and the rest of the neutron star are more efficient processes for the 
angular momentum transfer (Thompson et al. 2000). 



3. Discussion 

We have studied the propagation of shear- Alfven waves through the neutron-star crust to the 
exterior. Using a simplified cylindrical model, some interesting natures are found. The waves with 
high frequency, \uj/(mQ)\ ^ 1, always propagate, whereas those with low frequency propagate or 
non-propagate depending on spatial position. The transmission coefficient, which depends on the 
helicity state, is explicitly calculated. The ejection rate to the exterior increases with the frequency 
of the wave and the magnetic field strength. 

We now consider the behavior of modes in low frequency regime. The positive helicity modes 
£+ always propagate as far as c+ > 0, whereas the negative helicity modes £_ become evanescent 
near the surface for < er_ < (m + 2)Q. If corotating plasma covers the surface, the disturbances 
do not propagate in the atmosphere. On the other hand, if such plasma does not exist, i.e. static 
plasma or vacuum case, then the perturbations revive as outward waves in the exterior, although 
they suffer from some damping through the interior evanescent region. The exponential decay of 
the amplitude is not so large, for er_ ~ mil with large m, as shown in section 2.3. The disturbances 
therefore penetrate into the surface, and eventually escape to infinity. 

When the waves with the frequency < c_ < mil are ejected, they carry away negative 
angular momentum through the surface. In this case, the torque becomes positive and leads to 
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spin-up of the star by the outgoing waves. However, the magnitude is not enough to explain 
the glitches observed in AXPs(Kaspi et al. 2003). The geometrical factor may be important in 
actual situation. In our model, the spin axis agrees with the direction of wave propagation, and 
the direction of the perturbations is orthogonal to the spin axis. For magnetic field line mis-aligned 
with the rotation axis, the direction of the perturbations is perpendicular to the propagation, but 
not to the spin axis. In this case, 5B Z is induced in general, and the non-vanishing term 5B^,5B Z 
results in larger torque. See eq.(33) and the discussion. Compressional fast/slow magnetosonic 
waves are also induced at the same time. The calculations would be much more complicated. 

The increase of angular momentum as a result of radiation is very analogous to that in the 
gravitational radiation reaction instability, i.e. the CFS mechanism (Chandrasekhar 1970; Friedman 
& Schutz 1978). The generic criterion of the instability is o±lu < for unmagnetized rotating stars. 
Some fluid oscillation-modes satisfying this condition can grow in secular time scale as a reaction 
of gravitational radiation. The low frequency mode satisfies the instability criterion. Related 
with this, there is a suggestive work by Ho & Lai (2000). They considered the r-mode instability 
driven by Alfven wave emission as well as gravitational one. On the dimensional estimate, the 
transfer rate by Alfven wave is larger than that by gravitational wave, for highly magnetized star 
B > 3 x 10 12 (T/10ms)~ 3 G. They pointed out the Alfven wave-driven instability as an efficient 
process. In their treatment, however, the displacement current is neglected, and hence dynamical 
degree of the electromagnetic fields is eliminated in a sense. Such an approximation may be justified 
in the deep interior, but no longer near the surface. Our treatment is concentrated on the crust, 
but not the interior core, whose bulk motion is important to estimate gravitational radiation. Thus 
global calculation adequate for the emission of gravitational wave and Alfven wave is required to 
determine the efficient mechanism. 

As for a fundamental problem of the secular instability, the criterion should be clarified by 
taking into account of strong magnetic fields. For example, the canonical energy relevant to the 
instability is constructed only for unmagnetized stars (Friedman & Schutz 1978). In their treatment 
of the secular instability, the effect of gravitational emission is neglected in dynamical equation, 
but should be accounted in the energy equation governing the evolution in longer time scale. Our 
problem differs in this point. The degree of freedom of the electromagnetic fields is accounted in 
dynamical equation, and hence Alfven wave emission from surface is inevitably included. There 
is no mathematically rigorous proof of the secular stability criterion in magnetized stars, so that 
further study is necessary to conclude unstable growth of the modes coupled with electromagnetic 
fields. 

Finally, we will comment on some problems in applying to more realistic models relevant to 
magnetars. The rotation of observed magnetars is very slow Q ~ lrad s _1 . In order to satisfy the 
condition < cr_ < ml), a large azimuthal number, say, m ~ 10 3 is needed. If there are rapidly 
rotating stars, a moderate value of m may be sufficient. The other important factor is the magnetic 
configuration. Thompson et al. (2002) considered the magnetosphere threaded by large-scale electric 
current. The magnetic field structure is twisted, and quite different from our present model. The 
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perturbations of such magnetic fields are induced by bursts, and Alfven waves propagate. However, 
the detailed treatment seems to be complicated, and highly numerical calculations are needed. We 
expect that the behavior of short wavelength, say, much smaller than curvature of magnetic field, 
may be almost the same as given here. The result of the transmission rate calculated in section 
2.5 may be relevant. For the long wavelength mode, geometrical effects are important, and further 
studies are required. 
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